Structure formation in binary colloids 
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A theoretical study of the structure formation observed very recently [Phys. Rev. Lett. 90, 128303 
(2003)] in binary colloids is presented. In our model solely the dipole-dipole interaction of the 
particles is considered, electrohidrodynamic effects are excluded. Based on molecular dynamics 
simulations and analytic calculations we show that the total concentration of the particles, the 
relative concentration and the relative dipole moment of the components determine the structure 
of the colloid. At low concentrations the kinetic aggregation of particles results in fractal structures 
which show a crossover behavior when increasing the concentration. At high concentration various 
lattice structures are obtained in a good agreement with experiments. 

PACS numbers: 83.10.Pp, 82.70.Dd, 41.20.-q, 61.46.+W 
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Recently, the structure formation in a monolayer of 
colloidal particles under an external electric or magnetic 
field attracted considerable scientific interest. Investiga- 
tions have been performed with several different types 
of particles subjected to a constant or time dependent 
external helds. While in the absence of an external 
field chains, rings and fractal clusters have been ob- 
served with a temperature dependent fractal dimension 
0SSSEEHIESII3IIII in the presence of an ex- 
ternal field crystal structures emerg ed OOllSl 

It was 

demonstrated that by appropriate adjusting of the exter- 
nal field all the planar lattice structures can be produced 
in a magnetic monolayer ^1 > which is of high practical 
importance and also addresses essential problems of sta- 
tistical physics like the study of phase transitions between 
different morph olog ies and the phenomenon of melting in 
two dimensions |12lll3lll4| . 

Binary colloids, i.e. suspensions composed of two sorts 
of particles of the same size have been found very re- 
cently to produce novel structures when subjected to an 
ac electric field jl5j . Experiments performed with differ- 
ent frequencies and volume fractions of the components 
revealed a rich variety of self assembling structures in cer- 
tain frequency ranges: at frequencies where the effective 
particle interaction is repulsive a triangular lattice of par- 
ticles (Wigner crystal) was observed |l5|. At frequencies 
however where different types of particles attract each 
other, aggregation and crystallization occurred: at low 
concentrations strings and rings of alternating particles 
emerged while at higher concentrations various lattice 
structures were obtained ^5| . It was found that the mor- 
phology also depends on the relative concentration of the 
particles, i.e. if one type of particles has a sizeable excess 
flower-like clusters are formed with one type of particle in 
the middle surrounded by 5-6 particles of the other type. 
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In regions of 1 : 1 or 1 : 2 relative concentrations parti- 
cles crystallize into binary honeycomb or square-packed 
super-lattices. It has been shown that electrohydrody- 
namic (EHD) flow is alone not enough to understand the 
obtained structure formation but the dipolar interaction 
of the particles arising due to their induced dipole mo- 
ment might have a governing role to produce the observed 
structures |l5j . 

In this paper we present a theoretical study of the novel 
structure formation in binary colloids observed in Ref. 
[lol . Our goal is to understand the effect of dipolar forces 
on the structure formation, furthermore, to identify the 
governing parameters of the system and explore all pos- 
sible aggregate and crystalline morphologies determining 
the corresponding parameter regime of their occurrence. 

In the experiment the particles are settled down to the 
bottom plate of a container and an ac electric field is im- 
posed perpendicular to the bottom. Under the influence 
of the external field the particles attain a dipole moment 
parallel to the field, however, the magnitude and the ac- 
tual direction depends on the frequency of the field and 
can be different for particles made of different materials. 
In our model an ensemble of N particles of radius R is 
considered in a square box of side length L. Point-like 
dipole moments are assumed to be placed in the middle 
of the particles. For simplicity, the dipole moment of the 
particles is fixed during the time evolution of the system 
to be perpendicular to the plane of motion pointing either 
upward (particle type I) or downward (particle type II), 
which represents the two different material properties of 
the particles. This assumption implies that in the model 
only that frequency range of the experiment is considered 
where different particles attract each other. The parti- 
cles are considered to be suspended in an electromagnet- 
ically passive liquid so that EHD effects are completely 
excluded in the model. The liquid only exerts a friction 
force (Stokes force) on the particles which move under 
the action of dipole-dipole forces. When particles touch 
each other during their motion we introduce a repulsive 
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FIG. 1: (color) The aggregation of dipoles with fixed (j> r — 1, 
\jl t = 1 at different concentrations: a) (f> = 0.075, b) <f> = 0.125, 
c) (f> = 0.5, and d), = 0.70. For each simulation 1000 
particles were used and <f> was controlled by varying the side 
length L. 



force of the form of Hertz contact to prevent the overlap 
0, E| . After generating an initial configuration with ran- 
dom particle positions in the simulation box, and fixing 
the dipole moments according to the different material's 
properties, the time evolution of the system is followed 
by solving numerically the equation of motion of par- 
ticles for the two translational degrees of freedom with 
periodic boundary conditions (molecular dynamics simu- 
lation) . In the model the motion of particles is determin- 
istic, i.e. no stochastic forces are taken into account due 
to the high particle mass which hinders thermal motion 

15 L For details of our simulation techniques see Refs. 

3, li| • A model system of N particles consists of iVi and 
particles of dipole moment [i\ pointing upward (+), 
and dipole moment /i2 pointing downward (— ), respec- 
tively. The concentration of particles (j> is defined as the 
coverage <j> — NR 2 ir/L 2 . The partial concentration of the 
two components 4>i , 4>2 can be defined analogously whose 
ratio provides their relative concentration <p r — fa/fii- I n 
the model the magnitudes (Xi , [ii of dipole moments can 
be freely varied which captures up to some extent the 
effect of frequency tuning of the experiment. For sim- 
plicity, we fix [i\ and vary the ratio /i r = Two 
particles with parallel dipole moments exert a force onto 
each other which is isotropic (central), it always falls in 
the plane of motion parallel to the line connecting the 
two particles. The force has a 1/r 3 dependence on the 
separation r of the particles and it is repulsive for par- 
ticles of the same type and attractive for different ones. 
It is interesting to note that the system is rather simi- 
lar to an ensemble of charges where the interaction force 



FIG. 2: The number of particles n of clusters as a function of 
the radius of gyration R g . The dashed and continuous straight 
lines have slope 1 and 1.45, respectively. The crossover value 
n c was obtained as the crossing point of the straight lines 
fitted to the two regimes of different slopes. The inset shows 
n c {4>) where a straight line of slope 1.18 ± 0.05 was fitted to 
the data. 



decreases faster than the Coulomb force. 

Simulations have been carried out varying the three 
parameters of the model <p, <f> r and /i r in a broad range. 
Figure 2] shows simulation results obtained at different 
concentrations <fi fixing the relative concentration <f> r — 1 
and the relative dipole moment ji r = 1. It can be seen 
that at low concentrations the particles undergo a kinetic 
aggregation process and form chains and rings of alter- 
nating types of particles due to the attractive interaction 
of oppositely directed dipoles (Fig. ^1 and b) . It is inter- 
esting to note that the isotropic interparticle forces gives 
rise to anisotropic clusters since already a cluster of two 
particles behaves as an extended dipole oriented in the 
plane of motion which then results in a favorable direc- 
tion for the later aggregation steps. At low concentra- 
tions the chain is the dominating morphology, however, 
it can be shown analytically that for an even number of 
alternating particles the ring structure is energetically fa- 
vorable |l(| . (A ring can be observed in the upper part of 
Fig.m>.) The length of alternating chains is limited, when 
the chain length becomes comparable to the average dis- 
tance of chains, aggregation can occur not only at chain 
ends but also at internal particles resulting in branching 
and more compact structures. In order to characterize 
the structure of growing clusters we calculated the ra- 
dius of gyration R 2 g (n) = l/n(n — — ^i) 2 f° r 
each cluster during the time evolution of the system and 
averaged over clusters of the same size n. In Fig. [21 the 
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cluster size n is plotted as a function of the averaged 
gyration radius R g on a double logarithmic plot for sev- 
eral different concentrations. It can be observed in the 
figure that for each concentration n(R g ) is composed of 
two parts of power law functional form n ~ R™ : for small 
clusters n(R g ) can be well fitted by a power law of ex- 
ponent close to 1 indicating chainlike structures. At a 
certain chain length n c ((f>), however, a crossover occurs 
into a more compact structure characterized by a higher 
exponent of R g , a = 1.45 ± 0.05 was determined from 
simulations for [i r = 1 and <p r — 1. It follows from our 
argument that the crossover chain length n c has a power 
law dependence on the concentration, i.e. n c ~ (f>~" with 
an exponent (3 = 1.0. The numerical result is presented 
on the inset of Fig. where j3 = 1.18 ±0.05 was obtained 
in a reasonable agreement with the analytic prediction. 
Simulations also showed that the exponent a characteriz- 
ing the structure of growing aggregates does not depend 
on <j) until <j) < (jf ss 0.25. It is important to emphasize 
however that </v and /i r have a substantial effect on the 
aggregation process which will be explored in its entire 
complexity elsewhere |l6|. 

Increasing the concentration above <p* the morphol- 
ogy of clusters changes drastically; crystallites of square- 
packed structure are formed which then aggregate into 
large clusters. Fig. ^ shows simulation results obtained 
at <f> = 0.5 where a network of crystallites can be ob- 
served. At higher concentrations the entire system self 
assembles into a square-packed lattice of alternating par- 
ticles as in Fig. QJf. 

In order to understand the formation of crystal struc- 
tures in binary colloids and explore possible lattice mor- 
phologies beyond the square-packed one, we calculated 
the energy of basic particle configurations from which 
lattices can be built up. The Basic Configurations (BC) 
are composed of a particle of type I (with dipole moment 
/ii) surrounded by 2 (BCI), 3 {BCII), 4 (BCIII), 5 
(BCIV), or 6 (BCV) particles of type II (with dipole 
moment ^,2) placed at the energetically most favorable 
locations. In order to simply cover also cases where the 
different types of particles are interchanged, the energy 
of BCs was calculated as a function of \x r = ^2/1^1 fixing 
the value of /ii = 1 
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where (3 = n/(N — 1) and the particle number takes 
the values N = 3,4,5,6,7. The energy divided by the 
particle number E^(ii r )/N is presented in Fig. [3] as a 
function of fj, r for the possible values of N. It can be 
seen in the figure that all the BCs are stable (have a neg- 
ative energy) only in certain fj, r ranges between and a 
^} ax (N) since with increasing [i r the attraction exerted 
by the central particle is not enough to compensate the 
mutual repulsion of the surrounding ones. Such configu- 



FIG. 3: E N (^ r )/N from Eq. JTJ for N = 3,4,5,6,7. For 
clarity, the corresponding particle configurations are also pre- 
sented. 



rations where the central particle is surrounded by par- 
ticles of smaller dipole moment, i.e. fi r < 1 in Fig.^ are 
always stable except for BCV because y^ ax (7) w 0.8 
falls below 1. The upper bounds can be determined ex- 
actly from Eq. JTJ as fi™ ax (N) « 16,5.2,2.4, 1.3,0.8 was 
obtained for TV = 3, 4, 5, 6, 7, respectively. Repeating the 
above basic structures with alternating particles various 
structures can be built up, i.e. BCI results in chains or 
rings, BCII gives rise to a honeycomb lattice in which 
both types of particles have three neighbors of the other 
type, and BCIII leads to the square lattice. BCIV 
results only in flower-like structures, no lattice can be 
constructed since the plane cannot be covered by regular 
pentagons. Moreover, BCV forms the basis of a special 
type of honeycomb lattice, in which particles of the larger 
dipole moment have 6 neighbors of the smaller one but 
particles of the smaller moment have 3 neighbors of both 
types. The overall morphology attained by the entire par- 
ticle system at a given value of fj, r is not necessarily based 
on the energetically most favorable basic configuration. 
In spite of the isotropic particle-particle interaction the 
cluster-cluster interaction is highly anisotropic which can 
result in trapping the particles into local energy minima 
in configuration space, especially when thermal motion 
is hindered by the relatively large particle mass. Hence, 
the relative dipole moment fj, r , the total concentration (f> 
and the relative concentration <p r of components together 
determine the final structure. Based on the energy of ba- 
sic configurations Eq. (JTJ and Fig. [21 it is possible to 
determine analytically regions of </), (f> r , and \i T which are 
needed to obtain a certain lattice structure or aggregate 
morphology. The numerical results are summarized in 
Tableland Fig.0]provides an overview of possible struc- 



4 



FIG. 4: (color) Aggregate and crystalline morphologies ob- 
tained by simulations. In each system 1000 particles were 
used and concentration was controlled by changing the size 
of the simulation box. For visualization smaller pieces of the 
simulation box are cut out and magnified. 



tuxes obtained by simulations varying <f>, <fi r , \x r . 

If one of the components has a much larger concen- 
tration than the other one (for instance, <pi >> 4>2 ), 
due to the repulsive interaction of the identical particles 
a triangular lattice is formed with some binary islands. 
The limiting case when only one of the components is 
present is shown in Fig. 0Jt where a regular triangular 
lattice is observed. The simplest basic configuration BCI 
gives rise to chain and ring conformations which occurs at 
fi r = 1, 4> r = 1, in the concentration range (j) < (jf rs 0.25 
with a crossover into more compact morphologies as dis- 
cussed above. For an even number of particles the ring 
is energetically more favorable than the chain configura- 
tion. In spite of this, chains cannot close to form rings 
because the bending would require energy. Rings ap- 
pear when two nearby chains merge which occur most 
frequently at concentrations <f> w 0.125, see Fig.0J)). Sim- 
ilar chain and string formation of alternating particles 
has been reported in Ref. [l5j . 

At high concentrations lattice structures can be ob- 
tained, however, concentrations too close to the highest 
coverage (0 « 0.9) are disadvantagous because they pre- 
vent particle motion and result in freezing into local en- 
ergy minima. Hence, the best square lattice structures 
can be achieved at « 0.65 — 0.75. This type of lattice 
contains the same amount of both type of particles so 
that <f> r = 1 and fx r = 1 follows, see Fig. The param- 
eter regime providing square lattice agrees well with the 
experimental observation of Ref. 0] . The first type of 
honeycomb lattice (honeycomb I) where each particle has 
3 neighbors of the other type has not been observed ex- 
perimentally. However, our analytic calculations showed 
that energetically it can be favorable for the system in a 
certain parameter range. An example is presented in Fig. 
0Jf which was obtained at cf> = 0.5 and <p r = 1. For clarity, 



TABLE I: The parameter regimes determining final structure 
of the system. For explanation, see the text. (In the first 
row missing values of <f> r , \jl t indicate that solely one of the 
components is present.) 



Structure 





4>r 


u r 


triangular latt. 


0-0.9 






chain 


0-0.125 


1 


1 


ring 


0.125 


1 


1 


square lattice 


0.65 - 0.75 


1 


1 


honeycomb I 


0.45 - 0.50 


1 


« 2.5 


honeycomb II 


0.80 - 0.82 


2 


2.5- 


super lattice 


0.45 - 0.50 


« 2 - 2.5 


2.5-4.5 



some plaquettes of the lattice are highlighted in the fig- 
ure. The basis of this structure is BCII which is always 
stable for fi r < 1, i.e. when a larger dipole is surrounded 
by 3 smaller ones. However, to make the inverse config- 
uration more favorable, /i r has to be increased above the 
stability limit of BCII I (above ^ ax {5) w 2.4), other- 
wise, the system ends up in a square packed structure. 
fi r fss 2.5 proved to be an excellent choice numerically, see 
Fig. 0Ji If the magnitude of the two dipole moments are 
different (fi r ^ 1) concentration fluctuations easily lead 
to the formation of the second type of honeycomb lat- 
tice (honeycomb II), which has a hexagonal-closed-packed 
(HCP) structure, where particles of the larger dipole mo- 
ment have 6 neighbors of the other type and the ones 
with smaller dipole moment have 3 of both types. The 
highest portion of the system was found to crystallize 
into honeycomb II at « 0.80 — 0.82 with the ratio of 
4> r = 2 of the components and with high enough asym- 
metry of the magnitude of dipole moments fi r rs 2.5, 
which is needed to prevent the system to crystallize lo- 
cally into the square packed structure. The correspond- 
ing simulation results can be seen in F ig. life , which is in a 
nice agreement with the experiments 15]. If the relative 
dipole moment fj, r is even higher ( fi r > 2.5) but the con- 
centration and the relative concentration do not favor 
the emergence of honeycomb II, so-called super struc- 
tures can be observed in the colloid. These structures 
do not have ordered crystalline morphology, instead they 
are characterized by long straight binary chains which 
connect disordered or small crystalline island as in Fig. 

In summary, we proposed a simple model of a binary 
monolayer of dipolar particles to explore the effect of the 
induced dipole-dipole interaction in the intriguing struc- 
ture formation observed recently 0] . Varying the three 
parameters of the model, i.e. the total concentration 
of the particles </>, the relative concentration </> r and the 
relative dipole moment fi r of the components a rich vari- 
ety of structures were obtained in satisfactory agreement 
with the experimental findings 15]. The simplicity of 
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the model demonstrates that the main qualitative fea- 
tures of the structure formation observed are determined 
by the induced dipole-dipole interaction. Further theo- 
retical studies are in progress and further experiments 
are proposed. 
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